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<^ • We present new high resolution [R ~ 21, 000) near-infrared (A = 2/^m) spectroscopic 

. observations of a sample of Class I and flat-spectrum protostellar objects in the p 

Ophiuchi dark cloud. None of the five Class I spectra show CO v = - 2 absorption 



O \ features, consistent with high i^-band continuum veilings, 4 < < 20 and fast stellar 

\ rotation, assuming that the underlying protostellar photospheres are of late spectral 

' type, as is suggested by the low luminosities of most of these objects. Two of the 

flat-spectrum protostellar objects also show no absorption features and are likely to be 
^ ' highly veiled. The remaining two flat-spectrum sources show weak, broad absorptions 

which are consistent with an origin in quickly rotating {v sin i ~ 50 km s~^) late-type 
stellar photospheres which are also strongly veiled, ~ 3 — 4. These observations 
provide further evidence that: 1)-Class I sources are highly veiled at near-infrared 
wavelengths, confirming previous findings of lower resolution spectroscopic studies; and 
2)- flat-spectrum protostars rotate more rapidly than classical T Tauri stars (Class II 
sources) , supporting findings from a recent high resolution spectroscopic study of other 
flat-spectrum sources in this cloud. In addition our observations are consistent with 
the high rotation rates derived for two of the Class I protostellar objects in our sample 
from observations of variable hard X-ray emission obtained with the ASCA satellite. 
These observations suggest that certain Class I sources can rotate even more rapidly 
than flat-spectrum protostars, near breakup velocity. 



^Visiting Astronomer at the Infrared Telescope Facility which is operated by the University of Hawaii under 
contract to the National Aeronautics and Space Administration. 
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1. Introduction 

The physical natures, evolutionary states, and circumstellar disks of classical T Tauri stars 
are becoming better understood due to recent spectroscopic observations, high resolution imaging, 
and advances in the theory of pre-main-sequence (PMS) evolution. However, the natures of the 
central stars and inner circumstellar environments of protostars are still not very well known. 
This is primarily because they are so heavily extinguished that they are difficult to observe even 
with modern instruments and detectors. For example, it is not known whether the central stars 
in protostellar objects differ substantially in effective temperature, radius, or rotation properties 
from classical T Tauri stars (CTTSs). The presence of protostellar envelopes, significantly higher 
accretion rates, and more powerful outflows suggest that the photospheres of protostellar cores may 
indeed be physically different from CTTSs. Moreover, it is not known whether the physical natures 
of protostellar photospheres are consistent with the predictions of protostellar evolution and PMS 
stellar theory. 

Flat-spectrum and Class I young stellar objects (YSOs) are the best low-mass protostellar 
candidates for spectroscopic study because they have relatively well-developed central stars and 
are detectable at near-IR wavelengths. Several pioneering studies have recently been undertaken 
to begin investigating the physical natures of these objects. A few flat-spectrum protostars have 
been observed in low resolution spectroscopic surveys and thus far they appear to be characterized 
by late-type photospheres with high continuum veilings and sub-giant surface gravities (Luhman 
& Rieke 1999; Kenyon et al. 1998; Casah & Matthews 1992; Greene & Lada 1996, hereafter Paper 
I). Even fewer have been observed at high spectroscopic resolution, and these observations suggest 
that the flat-spectrum protostars rotate significantly faster than the more evolved CTTSs (Greene 
&; Lada 1997, hereafter Paper II). 

These investigations have provided some interesting clues to the natures of these objects, but 
more observations are clearly needed to make sense of these late-phase protostars. It would be most 
interesting to determine if the less evolved, more heavily embedded and veiled Class I protostars 
also show near-IR photospheric absorption lines when observed at high resolution with high signal- 
to-noise. Such observations could directly constrain the effective temperatures, gravities, veilings, 
and rotations of these objects, providing further evidence as to whether they are dominated by 
near-IR stellar, disk, or envelope emission and whether they are physically similar to flat-spectrum 
YSOs. More flat spectrum YSOs should also be observed with high-resolution near-IR spectroscopy 
to confirm that they are indeed late-type rapid rotators. In a sense the flat-spectrum YSOs provide 
a link to the well-known PMS stars, and this must be better developed so that they can in turn 
serve as a link to the less well-known Class I objects. 
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Therefore we have undertaken a new high-resolution, near-IR spectroscopic study of fiat- 
spectrum and Class I YSOs in the p Ophiuchi cloud core. We describe these new observations 
in §2 and present a rotation and veiling analysis of these data in §3. In §4 we discuss the likely 
natures of each of these objects and suggest further observational work. 

2. Observations and Data Reduction 

Near-IR spectra were acquired in 1997 May - June and 1998 July with the 3.0 m NASA Infrared 
Telescope Facility on Mauna Kea, Hawaii, using the CSHELL facility single-order cryogenic echelle 
spectrograph (Tokunaga et al. 1990; Greene et al. 1993). Spectra were acquired with a I'.'O (5 pixel) 
wide slit on the dates indicated in Table 1, providing a spectroscopic resolution R = X/SX = 21,000 
(14 km s^-*^). The spectrograph was fitted with a 256 x 256 pixel InSb detector array, and custom 
circular variable filters (CVFs) manufactured by Optical Coating Laboratories Incorporated were 
used for order sorting. These filters successfully eliminated the significant interference fringing 
normally produced in CSHELL and other echelle spectrographs which use CVFs for order sorting. 
The plate scale was 0'.'20 pixel" -"^ along the 30" long slit (oriented east - west on the sky), and all 
spectra were acquired at a central wavelength setting of 2.29353 fim corresponding to the v = 
- 2 CO band head. Each exposure had a spectral range A A ~ A/400 {Av ~ 700 km s"^). Total 
integration times for each YSO are given in Table 1. 

Data were acquired in pairs of exposures of up to 400 s duration each, with the telescope nodded 
10" east or west between exposures so that object spectra were acquired in all exposures. The BOV 
star HR 6165 (r Sco) was observed periodically for telluric corrections. The telescope was guided 
with the CSHELL internal CCD autoguider during exposures of these telluric correction stars, 
while the telescope tracking rates were adjusted for minimum drift while observing the optically 
invisible p Oph YSOs. Spectra of the internal CSHELL continuum lamp were taken for flat fields, 
and exposures of the internal CSHELL Ar and Kr lamps were used for wavelength calibrations. 

All data were reduced with IRAF. First, object and sky frames were differenced and then 
divided by flat fields. Next, bad pixels were fixed via interpolation, and spectra were extracted with 
the APALL task. Extracted spectra were typically 5 pixels (1") wide along the slit (spatial) direction 
at their half-intensity points. Spectra were wavelength calibrated using low-order fits to lines in the 
arc lamp exposures, and spectra at each slit position of each object were co-added. Instrumental 
and atmospheric features were removed by dividing wavelength-calibrated object spectra by spectra 
of early-type stars observed at similar airmass at each slit position. Final spectra were produced 
by combining the spectra of both slit positions for each object. 
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3. Data Analysis and Results 

3.1. Object Sample 

The object sample was selected from the Class I and fiat-spectrum YSOs observed at low 
spectral resolution (R ~ 500) in Paper I which were not subsequently observed at high spectral 
resolution in Paper ff . None of the newly observed sources (listed in Table 1) showed any absorption 
features in their low resolution iC-band spectra (Paper I), and they are also relatively bright, 7 mag 
< K < 10 mag. Table 1 shows that the brighter point sources (i.e. Elias 29, IRS 54) were observed 
with higher signal-to- noise ratios than the fainter ones (i.e. IRS 43, WL 6). This is useful for 
analyzing the veiling in these objects if the brightness differences among the sources are mostly 
due to different amount of IR excess emission from circumstellar regions. If this is true, then 
the brighter objects have greater near-IR veiling and greater signal-to-noise is required to detect 
their photospheric absorption lines. On the other hand, it may be possible that some of the bright 
protostars are featureless because they are of relatively early spectral type (G or earlier) and possess 
intrinsically weak K-band absorption lines (besides H and He). 

3.2. Veiling and Rotation Analysis 

The new fiat-spectrum and Class I YSO spectra are shown in Figures 1 and 2, respectively. 
None of the five Class I spectra show any evidence of CO absorption. Two of the flat-spectrum 
YSOs also show no evidence of any CO absorption (GSS 26 and YLW 13B), while the other two 
show evidence of weak, broad band heads and perhaps some overlapping rotation-vibration lines 
as well (Figure 1). 

We now analyze these spectra to constrain the veilings and physical natures of these sources. 
Our previous high resolution study (Paper II) showed that late-type Class II YSOs (pre-main- 
sequence stars) rotate slowly, (vsini) < 20 km s~^, while fiat-spectrum ones rotate quickly, {vsini) > 
20 km s~^. This study also showed that in addition to broadening the band head, high rotation 
also decreased the maximum absorption depth of the band head and adjacent individual rotation- 
vibration lines. Thus rotation as well as veiling can reduce the detectability of CO absorption in 
finite signal-to-noise spectra. 

We now explore the limits of i^-band veiling and rotation in our new protostar sample by 
comparing their spectra to those of flat-spectrum and Class II YSOs to which we have artiflcially 
added continuum veiling. We chose VSSG 25 to be representative of a slowly rotating Class II 
YSO {v sin z = 5 km s~^) and VSSG 17 to be representative of a quickly rotating {v sin i = 47 
km s~^) flat-spectrum source (see Paper II). The spectral types and X-band veilings have been 
measured for both of these objects. VSSG 25 has a spectral type MOIV/V with = 0.25 (Luhman 
& Rieke 1999) while VSSG 17 is MOIV/V with rfc = 0.9 (Luhman & Riekc; Paper I). The i^-band 
veiling is defined as = Fkcx/Fk* where Fkcx is the K band excess fiux, and Fk* is the K band 
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stellar flux. We added veiling (a constant positive offset) to the template spectrum of VSSG 17 

(taken from Paper II) and degraded its signal-to-noise so that its CO absorption equivalent width 
and detectability were weakened. Wc then used a series of these VSSG 17 templates with various 
veilings and signal-to-noise ratios to estimate the veilings of our observed sources. 

We estimated the veiling for sources which show CO band head absorptions (IRS 51 and IRS 
63) by matching their spectra to veiled templates with identical signal-to-noise and similar features 
(CO band head depth and slope). The band head profiles of these objects also matched the shape of 
VSSG 17 and matched those of observed slowly rotating late-type stars which had been artificially 
broadened with a stellar rotation profile of v sin i = 50 km s""^ (see Paper II). Thus it is likely 
that their CO absorptions arise in rapidly rotating stellar photospheres. However, the absorption 
features of these highly veiled sources are very weak, and this limits the derived rotation velocities 
to uncertainties of approximately 40% and the derived veilings to uncertainties of at least 20%. 

The minimum likely veilings of the Class I and flat-spectrum sources without detectable CO 
absorptions were estimated by assuming that their CO absorptions were intrinsically similar to 
VSSG 17 but had increased continuum veiling. We derived an analytical relation for veiling based 
on the principle that an object's CO band head is undetectable when its maximum absorption depth 
is less than 3.0 times the RMS noise over 1 resolution element (5 pixels) in the spectrum. This 
is a robust minimum likely veiling criterion because less veiling would ensure that the band head 
would be definitely detected, and considering more resolution elements would increase the amount 
of veiling derived. The resultant derived minimum likely values matched those of templates 
which were veiled to the point where their CO absorptions just disappeared visually. 

Possible differences between the actual spectral types of these sources and the MO template 
also cause uncertainties in the derived veilings of up to about 50%, but most deviations are expected 
to be smaller than this. All observed Class I sources except Elias 29 are relatively low luminosity, 
2.4 Lq < Lboi < 13 L0 (Wilking, Lada, k Young 1989, hereafter WLY). All observed fiat-spectrum 
YSOs have 1 Lq < Lboi < 3 Lq (WLY; Greene et al. 1994). These luminosities are consistent with 
the sources being low-mass PMS YSOs (M < 1 Mq) if the Class I objects are powered by accretion, 
an assumption which is consistent with their high veilings (see also §4.2). In using VSSG 17 as a 
spectral type template to measure veilings, we have implicitly assumed that all sources have spectral 
types near MOIV/V, a typical value for T Tauri stars. A young star of this spectral type has a mass 
of approximately 0.4 M© if on the birthline of the H-R Diagram (sec Stabler 1988). If an observed 
YSO is really 1 Mq, then it would have a spectral type of K3-4IV/V (see D'Antona & Mazzitelli 
1997; Stabler 1988) and thus an intrinsic stellar CO equivalent width of only about 0.6 times that of 
an MOIV/V star (see Paper I). This would mean that its actual veiling is rj^ = 0.6rfe — 0.4 (see §4.2 
of Paper I) where is the veiling estimate made assuming an MOIV/V spectral type. Likewise, if 
an observed object were really a lower mass star near the brown dwarf limit, then its true birthline 
spectral type would be near M6IV/V and its intrinsic CO absorption would be about 40% greater 
than that of the MOIV/V template. Thus the true veilings of observed sources may differ from 
our derived ones by as much as 50%. However, it is unlikely that many of the observed objects 
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have such large deviations (in either direction) because their luminosities are consistent with most 

having masses of approximately 0.5 Mq, similar to that expected for our MOIV/V template. We 
have used other published information on the observed objects when available to estimate their 
intrinsic spectral types before deriving veiling estimates (see §4.1 and 4.2). 

Finally, we studied how rotation decreases the apparent CO band head absorption depths. 
This effect can be separated from veiling by measuring the slope of the CO band head, but this 
is not possible for sources which do not show this absorption feature. The majority of sources in 
our sample do not show this feature, so we do not know their rotation velocities. Our earlier study 
(Paper II) showed that fiat-spectrum sources rotate significantly faster than Class II YSOs, so our 
current sample of more highly embedded flat-spectrum and Class I YSOs maj^ rotate even faster 
still (perhaps v sin i > 50 km s~^). Indeed, recent X-ray observations with the ASCA satellite 
indicate that one Class I source in our sample (IRS 43) is rotating at least this rapidly (Montmerle 
et al. 2000). We estimate that the equatorial rotational breakup velocities of these young stars 
(M ~ O.5M0 and R ~ 2Rq) are approximately v 2± 180 km s"^, or v sin z ~ 150 km s"^ for a 
mean inclination i = 57deg. Next, we studied how rotation decreases the maximum CO band head 
absorption in YSOs by artificially rotating our WL 5 template by convolving its observed spectrum 
with limb-darkened stellar broadening profiles for 25 km s~^ < v sin i < 175 km s~^. See Paper II 
for more details and examples of artificially rotated spectra. 

These experiments showed that the maximum CO absorption depth of the VSSG 25 template 
with V sin i = 150 km s~^ is a factor of 1.37 weaker than the one with v sin i = 50 km s~^. 
Therefore if VSSG 17 were rotating at breakup, its maximum CO absorption depth would be only 
73% as deep as now seen in its spectrum. Thus if rotating near breakup, the featureless objects 
in our sample would have somewhat lower continuum veilings than those calculated based on the 
observed VSSG 17 spectrum. We have calculated these reduced values, and in Table 2 we present 
all of these estimated veilings for the observed sources, using our best estimates of their intrinsic 
spectral types. 

4. Discussion and Conclusions 

We now discuss how the results of this veiling / rotation analysis and pre-existing data constrain 
the possible physical natures of these sources. 

4.1. Flat-Spectrum Objects 

The flat-spectrum YSOs IRS 63 and IRS 51 were both found to have broad, weak CO absorp- 
tions which matched those expected for late-type stellar photospheres rotating at v sin z ~ 50 km 
s~^. The weak CO absorptions of these two YSOs are consistent with their not being detected in 
our initial low-resolution survey (Paper I) . We estimate the continuum veiling of IRS 63 to be ~ 4 
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provided that it is a PMS YSO near MO spectral type. Luhman k, Rieke find that the spectral 
type of IRS 51 is G5 - K7, earlier than our MO template VSSG 17. Thus their derived veiling = 
1 - 3 is lower than ours because a G5 - K7 PMS star has less intrinsic CO absorption than a MO 
one (see §3.2). It is likely that IRS 51 is indeed an embedded low-mass YSO because its bolometric 
luminosity is only 1.4 Lq (WLY). The birthline mass for this luminosity is approximately 0.5 Mq, 
corresponding to a spectral type of K5-7 and a true veiling of — 3. Thus it is likely that both 
IRS 63 and IRS 51 are similar to the quickly rotating flat-spectrum YSOs which we analyzed in 
Paper II, but these new objects have even greater veiling (i.e. = 3 - 4 versus V}^ 1 for the 
Paper II YSOs). 

Luhman & Rieke found GSS 26 to have variable veiling, = 0.75 and = 4 at epochs of 
1994 July and 1996 May, respectively. Our spectrum of GSS 26 (in Figure 1) was taken in 1997 
May, one year after the latest Luhman k. Rieke spectrum. We estimate that > 11 when our 
spectrum was acquired, and our assumption of a MO spectral type is consistent with the Luhman & 
Rieke determination of K5 - M2. This rapid increase in veiling - a factor of 2 each year - is perhaps 
suggestive of a similarly rapid increase in accretion. Luhman &; Rieke also note that this source 
increased in brightness by K ~ 1.2 mag between epochs. Our spectra are not photometrically 
calibrated, but comparisons with other objects support that this source was at least as bright as 
when observed by Luhman & Rieke the previous year. This is one of the YSOs Luhman &; Rieke 
observed whose HI Br 7 emission line flux increased as its veiling increased, implying that the excess 
continuum emission is associated with a circumstellar accretion disk if the Br 7 emission arises from 
disk accretion. The rapid variability of this object's veiling also suggests that its if-band veiling 
is produced by accretion from inner disk distances (several AU) and not from an outer disk or an 
outer circumstellar envelope. 

The high values of these objects also constrain the physical origins of their veilings. In Paper 
I we showed that veilings > 1 cannot be produced by a simple optically thick, geometrically thin 
reprocessing disk around a low-mass PMS star. Consequently we argued that these high veilings 
are most likely produced by either actively accreting circumstellar disks or circumstellar envelopes 
associated with these objects. Furthermore we found that veilings in the range measured for IRS 
51 and IRS 63, = 3 - 4, can be caused by luminous accretion disks {LfHsk/L* < 3). Veilings in 
the range observed for GSS 26, ~ 5 — 10, could be explained by extremely luminous accretion 
disks [Ldisk/L^: > 3). In either case, these accretion disks would have to have relatively large 
central holes to avoid producing strong CO absorption-line systems in the disk photosphere itself. 
However, because there is no obvious physical mechanism for producing central holes of the needed 
size, Calvet, Hartmann, & Strom (1997) suggested that the veiling flux must originate in some other 
circumstellar structure such as the inner regions of the protostcllar envelope. On the other hand, 
this seems to be inconsistent with the observation by Luhman & Rieke that the i^-band excesses 
of flat-spectrum YSOs arc correlated with their HI Br 7 line fluxes which in turn suggests that 
the veiling flux should originate in the disk. More detailed knowledge of the conditions required to 
produce CO absorption line systems in an accretion disk may be needed to resolve this issue. 
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The flat-spectrum source YLW 13B was found to have H Br 7 absorption by Luhman & Rieke, 
who estimate its spectral type to be earher than KO. However, they also find it to be significantly 
veiled with > 1, so it is possibly an intermediate mass PMS cloud population member. Our 
non-detection of CO absorption does not constrain this source further. 

4.2. Class I Objects 

We do not detect CO absorptions in any of the Class I YSOs which we observed (Figure 
2), confirming earlier low resolution spectroscopic observations that found all these objects to be 
featureless and highly veiled (Paper I; Luhman & Rieke). Consequently their spectral types are 
unknown, however all of these Class I sources exhibit HI Br 7 emission (Paper I; Luhman & Rieke). 

Analysis of our new data constrains the natures of these objects. Table 2 shows that they all 
have large veilings, > 4, if they are late-type (~ MO) stars. GSS 30, IRS 43, and WL 6, all 
have estimated minimum veilings of — 5 — 8, overlapping with the flat-spectrum sample. These 
sources all have bolometric luminosities Lboi < 13 Lq (WLY). This is consistent with their being 
low-mass (M < 1 Mq) protostars accreting matter at rates M ~ 5 x 10~® Mq yr~^, the value 
expected for the T ~ 20 K gas temperatures in the p Oph cloud (see Adams, Lada, & Shu 1987, 
hereafter ALS). Veilings in the observed range, > 4 — 10, are also predicted by theoretical models 
of Class I circumstellar envelopes (Paper I; Calvet et al.). 

Montmerle et al. note that the maximum possible mass of IRS 43 (also known as YLW 15) 
is 2.2 Mq which is derived from PMS models given its bolometric luminosity (L ~ 10 L©) and 
assuming it is on the birthline. Likewise, they calculate that the maximum likely mass of WL6 is 
approximately 0.4 Mq; this increases somewhat if the Lboi = 2.4 Lq of WLY is adopted. These 
maximum calculated masses assume that essentially all luminosity is due to photospheric thermal 
radiation (L = A^irR^aT'^) and essentially none is due to accretion [L = GMM/R). These sources 
are discussed further in §4.3. 

We estimate that Elias 29 and IRS 54 have very high veilings if they are late- type low-mass 
stars, rfe > 14 — 34. The bolometric luminosity of IRS 54 is estimated to be only 12 Lq (WLY), 
also consistent with this object being a low-mass protostar which is accreting its envelope at the 
rate prescribed by the p Oph cloud's gas temperature. However, its continuum veiling must be 
Tfc > 14 — 20 if it has spectral type MO and is rotating rapidly. This is about a factor of 2 higher 
than the model predictions of Calvet et al., but those calculations were done for a hypothetical p 
Oph Class I YSO with L = 5 Lq. The model may predict greater veiling for IRS 54 if its higher 
luminosity is taken into consideration. It is also possible that this source may be a somewhat earlier 
type protostar which is less veiled. 

Elias 29 has the highest derived veiling of the sample, rj^ > 25 — 34, assuming an intrinsic MO 
photosphere and a high rotation rate. WLY estimate its luminosity to be Lboi = 48 Lq, and ALS 
have modeled it as a 1 Mq protostar which is accreting its circumstellar envelope. Such a star 
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would have a spectral type of K3-4 if on the birthhne (see D'Antona & MazzitelH 1997; Stahler 
1988), with an intrinsic CO absorption approximately 60% as strong as that of an MO star (see 
§3.2). Thus wc revise our estimate of the likely veiling of this YSO to > 15 — 20 if it is indeed a 
1 Mq protostar. In Paper I we analyzed the ALS model for Elias 29 and showed that the predicted 
emission from the inner protostellar envelope of this source would produce a veiling of 20, 
assuming that its disk luminosity is 0.75 Lboi- Thus our new measurement is consistent with our 
earlier prediction based on the ALS model. 

The high luminosity of Elias 29 also allows for it being a more massive, earlier spectral type 
YSO that has higher stellar luminosity and less accretion luminosity than assumed by the ALS 
model. However, it is unlikely to be very different because the observed near-to-far IR energy 
distribution and the 10 ^m silicate absorption of Elias 29 are fit well by the ALS model, and there 
are no clues which indicate that Elias 29 is an early-type object. For example, WL 16, which is likely 
an early A type star (Biscaya Holzbach et al. 2000), has mid-IR aromatic hydrocarbon emission 
features which indicate a UV radiation field (Hanner, Tokunaga, & Geballe 1992). However, Elias 
29 shows no evidence for IR hydrocarbon emission (Hanner, Brooke, Sz Tokunaga 1995; Boogert 
1999) and thus no evidence for a UV radiation field. 

4.3. Protostellar Rotation 

Our observations, specifically the broad band head shapes of IRS 51 and IRS 63, strengthen 
the earlier findings of Paper II which suggested that flat-spectrum protostars rotate more rapidly 
than Class II sources (CTTSs). Since flat-spectrum protostars are believed to be evolutionary 
precursors of CTTSs, then this finding may indicate that certain physical conditions characteristic 
of protostellar evolution (e.g., high accretion rates) may result in their higher rotation rates. It is 
therefore interesting to ask whether the less evolved Class I sources might rotate even more rapidly 
than the flat spectrum sources. Indeed, some models of protostar development predict that such 
objects should be rotating near breakup (Shu 1991). 

Recently, strong hard X-ray flares have been observed with the ASCA satellite from two of 
the Class I protostars in our sample - IRS 43 and WL 6 (Montmerle et al.). These remarkable 
observations reveal relatively rapid periodicities in the X-ray emission from these sources, enabling 
the derivation of their photometric rotation rates. Montmerle et al. find that WL 6 is rotating with 
a period of about 3 days {v sin z 2± 40 km s~^ for a 0.5 Mq star on the birthline), comparable to 
the rotation rates of the flat-spectrum sources observed here and in Paper II. This source also has 
a weak outflow (Sekimoto et al. 1997), undetected millimeter emission from its envelope (Andre & 
Montmerle 1994), and an IR energy distribution which can be modeled as a highly extinguished 
flat spectrum YSO (Montmerle et al.). All of these properties indicate that WL 6 may indeed be 
very similar to the fiat spectrum YSOs for which we have detected absorption lines and have found 
to be rotating more rapidly than CTTSs (Class II YSOs). 
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Montmerle et al. also argue that the central star of IRS 43 has a 20 h rotation period, the 
observed period of its X-ray variability. This requires that its mass be greater than or equal to 1.8 
Mq in order for it to be rotating below breakup velocity if it is on the birthline (Montmerle et al.). 
Thus the mass of IRS 43 is constrained to lie in the range 1.8 - 2.2 Mq by its X-ray emission and 
bolometric luminosity (see §4.2) if its rotation period is indeed equal to its X-ray variability period 
of 20 h. IRS 43 does have a more steep mid-IR energy distribution (clearly Class I) than WL 6, and 
it also has spatially-resolved (r ~ 3000 AU) millimeter emission with a derived envelope mass of 
approximately 0.1 M© (Andre &; Montmerle). Therefore it is likely to be in an earlier evolutionary 
state than WL 6. 

The slow rotation velocities of CTTSs have been explained by angular momentum regulation 
of these stars by magnetic coupling to their disks. Edwards et al. (1993) found that late-type T 
Tauri stars (TTSs) with large H — K IR excesses (CTTSs) had slow rotation periods, P > 4 d. 
They also found that TTS with small H — K IR excesses had a broad range of periods, including a 
significant number with P < 4 days. Edwards et al. interpreted this correlation to arise because the 
magnetic fields of the CTTSs were were coupled to their disks, providing stellar angular momentum 
regulation and therefore long stellar rotation periods. The low-excess TTSs had already dissipated 
their disks and so were not subject to this regulation mechanism. More recent studies of larger 
TTS samples have both disputed that the correlation between IR excess and rotation period exists 
(Stassun et al. 1999) and have provided evidence that it exists but is weak (Herbst et al. 2000). 

Our studies of YSO rotation (this paper and Paper II) have shown that flat-spectrum protostars 
rotate significantly more rapidly than Class II YSOs or CTTSs, suggesting that rotation velocities 
decrease as stars evolve past the protostellar state. This scenario has been bolstered and expanded 
further by the recent X-ray results of Montmerle et al.. Taken together (and along with the many 
rotation studies of optically visible CTTSs), these works suggest that heavily embedded protostars 
(Class I) rotate very rapidly, in some cases near breakup velocity, while less embedded ones (flat- 
spectrum YSOs) rotate somewhat less rapidly, at about 1/3 breakup velocity {v sin i ~ 50 km s~^), 
and Class II YSOs / CTTSs rotate slowly, v sini < 20 km s~^. This finding would be strengthened 
considerably by further cross-checking of observational techniques; the X-ray protostars should be 
observed at higher signal-to-noisc in the ncar-IR to search for rotationally broadened lines, while 
the flat spectrum and Class II YSOs with IR-derived rotation velocities should be observed for 
periodic X-ray variability. 

If the angular momenta of low-mass YSOs are indeed regulated by star— disk coupling, then 
the fact that flat-spectrum (and at least one Class I) YSOs rotate significantly more rapidly than 
CTTS implies that either the flat-spectrum / protostellar YSOs are coupled to faster rotating 
disk regions than CTTSs, or else that stars and disks do not become rotationally locked until the 
CTTS evolutionary phase. In the first case, flat-spectrum and Class I YSOs may couple to their 
disks at smaller radii (and hence have higher rotation velocities) because their accretion rates are 
much higher than CTTS. The veilings and luminosities of Class I and flat-spectrum YSOs are 
considerably higher (by about an order of magnitude) than CTTS, supporting the notion that they 
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have higher accretion rates also. In support of the second case, Montmerle et al. have posited that 
protostars arc not initially magnetically coupled to their disks but rather spin-down and become 
coupled over a magnetic braking time on the order of 10^ yr which is nearly linearly proportional 
to stellar mass. This is comparable to the lifetime of the Class I and flat spectrum phases, so this 
would account for the higher rotation velocities of flat-spectrum and Class I protostars. This latter 
magnetic braking scenario of velocity evolution from Class I to fiat-spectrum to CTTS YSOs may 
be somewhat complicated by mass effects; Montmerle et al. predict that at the same age more 
massive protostars will rotate more quickly than less massive ones. Protostellar masses must be 
measured much more accurately before this effect can be verified, however. 

Obtaining new high resolution, high signal-to-noise spectra over the entire 1.5 — 2.4 jum region 
will likely be the best method for obtaining more definitive information on the masses (spectral 
types) and rotational characteristics of Class I protostars. This wide spectral range is required in 
order to be sensitive to a wide range of spectral types. Even intermediate-to-high mass stars with 
HI line emission may show H-hand HI Br absorption lines which may strongly constrain spectral 
types and masses (e.g., Biscaya Holzbach et al.), while -RT-band data are required to determine the 
properties of very red late-type YSOs. Detecting and resolving near-IR lines in WL 6 and IRS 43 
would allow determination of their masses and photosphcric rotation rates, providing a good test 
of the emerging scenario of protostellar rotational evolution. 

We thank the referee John Lacy and also Pat Cassen for providing comments which improved 
this paper. We also thank W. Golisch, D. Griep, and C. Kaminski for assistance with the ob- 
servations. We acknowledge the National Science Foundation for funding grant AST-9420506 to 
develop the fringe-free CVFs used to acquire these data with CSHELL. TPG acknowledges a grant 
from the NASA Ames Research Center Director's Discretionary Fund. All data were reduced with 
IRAF, which is distributed by the National Optical Astronomy Observatories, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under contract to the National 
Science Foundation. 
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Table 1. IRTF CSHELL 2.2935 /xm CO Observations 



Source 


a(1950) 


5(1950) 


UT Date 


Int. Time 


S/N^ 


Features 




hh mm ss.s 


o / // 




(minutes) 






p Oph Class I YSOs: 












GSS 30 


16 23 20.0 


-24 16 22 


1999 Jul 07 


74.0 


25 


none 


Elias 29 


16 24 07.8 


-24 30 35 


1997 Jul 07 


46.0 


135 


none 


WL 6 


16 24 20.0 


-24 23 11 


1998 Jul 04 


102 


25 


none 


IRS 43 


16 24 25.1 


-24 34 10 


1997 May 05 


56.7 


35 


none 


IRS 54 


16 24 50.0 


-24 25 05 


1997 May 05 


63.3 


80 


none 


p Oph Flat- 


■Spectrum YSOs: 










GSS 26 


16 23 09.0 


-24 14 11 


1997 May 5 


30.0 


65 


none 


YLW 13B 


16 24 19.3 


-24 35 03 


1997 May 5 


40.0 


65 


none 


IRS 51 


16 24 37.6 


-24 36 35 


1997 Jul 4 


60.0 


50 


weak broad abs 


IRS 63 


16 28 34.3 


-23 55 05 


1997 Jul 4 


43.3 


40 


weak broad abs 



^Approximate signal-to-noise ratios were determined in relatively featureless regions of 
spectra in regions of good atmospheric transmission. 
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Table 2. Derived Veiling Estimates 



Source r^* ^'fe.break'^ Comment 



p Oph Class I YSOs: 



GSS 30 


5.5 


3.8 


lower limits 


Elias 29 


20^= 


15^= 


lower limits 


WL 6 


5.5 


3.8 


lower limits 


IRS 43 


8.1 


5.7 


lower limits 


IRS 54 


20 


14 


lower limits 


p Oph Flat- 


Spectrum 


YSOs: 




GSS 26 


16 


11 


lower limits 



YLW 13B • • • • • • 

IRS 51 3*^ measured value 

IRS 63 4 measured value 



'^Minimum likely /f-band veiling assuming intrinsic CO depths and rotation 
velocities similar to VSSG 17. See §3.2 for uncertainties. 

'^As for note (a) but computed for sources rotating near breakup velocity. 

^Thc minimum veilings shown were computed for a 1 Mq object, 
with intrinsic CO absorptions 40% weaker than VSSG 17. (see §4.2 

■^Wc do not report our derived veilings since this is clearly not a late-type 
object (Luhman & Rieke 1999; see §4.1). 

^The Luhman & Rieke (1999) spectrum and Lboi = 1-4 of IRS 51 con- 
strain its spectral type to K5 - K7, which we have assumed to derive its veiling 
(see §4.1). 
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p Ophiuchi Flat Spectrum YSOs 



X 

> 




GSS 26 



IRS 51 



_ IRS 63 



2.291 2.292 2.293 2.294 2.295 

Observed Wavelength (|im) 



■ - YLW 13B 



2.296 



Fig. 1. — Spectra of flat-spectrum p Oph protostars in the v = - 2 CO band head region. 
The approximate location of the 2.2935 /xm band head is indicated (BH), and the approximate 
continuum level of the region red-ward of the band head is shown by a dotted line for each spectrum. 
Each spectrum has been scaled to a mean value of 1. YLW 13B has a zero point of 0, but the other 
spectra are offset by a value of 0.2 relative to each other. The sources IRS 51 and IRS 63 show 
clear evidence for CO band head absorptions, while GSS 26 and YLW 13B do not. The spectral 
resolution element is approximately 1.1 x 10~^/xm, so the band heads of IRS 51 and IRS 63 are 
clearly broadened by rotation. 
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p Ophiuchi Class I YSOs 



1.^ 



1.6 



I 

> 



1.4 



C5 1-2 



1.0 




- - EUas 29 




GSS 30 



IRS 43 



IRS 54 



WL6 



2.291 



2.292 2.293 2.294 2.295 

Observed Wavelength (|am) 



2.296 



Fig. 2. — Spectra of Class I p Oph protostars in the v = - 2 CO band head region. The ap- 
proximate location of the 2.2935 /xm band head is indicated (BH), and the approximate continuum 
level of the region red-ward of the band head is shown by a dotted line for each spectrum. Each 
spectrum has been scaled to a mean value of 1. WL 6 has a zero point of 0, but the other spectra 
are offset by a value of 0.2 relative to each other. No Class I spectra show evidence of any CO 
absorption. 



